pub fn new_birthday_probability(n: u32) -> f64 {
    let mut res = 1_f64;
    for i in (365 - n + 1)..=365 {
        res *= i as f64 / 365_f64
    }
    1_f64 - res
}
